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We show that weU-posed, conformaUy-decomposed formulations of the 3+1 Einstein equations can 
be obtained by densitizing the lapse and by combining the constraints with the evolution equations. 
We compute the characteristics structure and verify the constraint propagation of these new well- 
posed formulations. In these formulations, the trace of the extrinsic curvature and the determinant 
of the 3-metric are singled out from the rest of the dynamical variables, but are evolved as part of 
the well-posed evolution system. The only free functions are the lapse density and the shift vector. 
We find that there is a 3-parameter freedom in formulating these equations in a well-posed manner, 
and that part of the parameter space found consists of formulations with causal characteristics, 
namely, characteristics that lie only within the lightcone. In particular there is a 1-parameter family 
of systems whose characteristics are either normal to the slicing or lie along the lightcone of the 
evolving metric. 

PACS 



I. INTRODUCTION 



, Analytical work in recent years has produced a number of systems of evolution equations which are equivalent to 
the Einstein equations at the constraint manifold, and which have a well posed initial value formulation |^^. 
^) . What motivates interest in this type of result is a general understanding (see for instance [||) that explicit well- 
posedness would be relevant in implementing consistent and stable numerical algorithms to integrate blackhole space- 
times. 

The well-posed schemes for which a numerical code has been implemented appear not to exibit significant improvc- 
■ ments over other methods, there being several factors relevant to numerical implementation which play a significant 
O ■ role. What is puzzling, however, is that, on the other hand, there have been numerical simulations with apparently 
^ '■ better behavior, but which are based on systems which do not seem to have the well-posed character. One preponder- 
. ant feature of these numerically more robust schemes is that they are built on a decomposition of the intrinsic metric 
into a metric of unit determinant and the determinant itself, and of the extrinsic curvature into trace and trace-free 
part. With slight variations, this way of evolving the 3-1-1 Einstein equations has been considered by [pl jlC | . Quite 



, , ■ recently, this form has been shown to possess striking computational advantages over the standard form We 
rN \ refer to this general scheme as a conformally-decomposed formulation of the 3-1-1 Einstein equations. 
^ . It is difficult to explain the success of these systems as opposed to the well-posed evolution schemes, or even to the 
standard (ADM) evolution schemes. The relative sizes of the fields of well-posed systems are roughly the same for 
different spectral frequencies in a Fourier representation, which helps explain the stability of the system via numerical 
analysis. However, in the conformally-decomposed systems this does not happen in general (for standard norms), as it 
does not happen for the standard ADM system, thus making it more difficult to justify their relative better behavior. 

We can speculate on two features that can possibly bear relevance to well behaved numerical evolution. One feature 
is that good behavior in evolution is related to constraint violations. If the system preserves more accurately the 
constraints, then the evolution remains closer to the constraint submanifold, which contains the physical solutions. 
Outside this submanifold the solutions are unphysical; thus, there is no compelling reason to rule out fast growths 
for seemingly tame initial data for unphysical solutions. Thus, we suggest that controlling the constraint violations 
may lead to well-behaved numerical evolution. In this respect, it has been shown |l2), that (at least in the linearized 
case) there are well-posed modifications of the Einstein equations outside the constraint submanifold which make that 
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submanifold an attractor, thus improving the chances of building numerical codes with better behaved constraint 
propagation. Another possible cause of concern for generating numerical instabilities is the nature of the boundary 
conditions which are ussually imposed. The initial-boundary value problem for the Einstein equations has not been 
studied for the systems used in numerical simulations and where instabilities have been found (see nevertheless [p^ for 
a complete theory of boundary values for Einstein's equations in conformal frame variables, and p4[ for a linearized 
study), thus the set of boundary conditions for which the constraint equations are satisfied is in general not known. 
In dealing with this problem, establishing well-posedness for the Cauchy problem is a necessary first step. 

The other feature which could give rise to numerical instabilities is the relative sizes of the "longitudinal" and the 
"radiative" modes in general relativity. In all non-trivial asymptotically flat solutions (either vacuum or with matter 
satisfying the appropiate energy conditions) the positivity of the mass implies the existence of longitudinal modes, 
and there are many astrophysically relevant cases where there is an approximate local notion of longitudinal versus 
transverse modes, and where the former are several orders of magnitude bigger that the latter. If they are not properly 
separated in the numerical algorithms, the errors caused by finite differencing might be of the order of the "radiative" 
modes, and bad behavior can be expected. The separation of the conformal freedom in the more successful codes can 
perhaps be thought of as a way of dealing with this issue, or at least isolating it. 

In this work, we focus on this latter aspect. A technique for taking advantage of the conformal factor to partially 
decouple the "longitudinal" and "transversal" modes was used to obtain results on the Newtonian limit of general 
relativity [^j. In that case the conformal field was fixed via an elliptic equation, decoupling in this way the more 
prominent Newtonian potential to first order from the radiative degrees of freedom. Further studies on this problem 
would be critical to obtain realistic simulations of most astrophysically relevant problems. 

Here we construct 3-parameter families of first-order well-posed systems wich share some of the properties of the 
more successful systems, such as the conformal decomposition of the fundamental fields, in the hope that their study 
would help understand what is causing them to behave better than others. In Section || we apply techniques similar 
to those we used in [^|j7j in order to obtain versions of the 3-1-1 equations that are conformally-decomposed but which 
are well posed. Additionally, we calculate the structure of characteristics and show that for a open region in parameter 
space the resulting equations are metric-causal, namely they have all propagation cones inside or coincident with the 
light cone. There is even a one parameter subfamily for which propagation is either along the light cone or normal to 
the slices. 

Furthermore we show that the constraints are propagated by these well-posed evolution equations. As oposed 
to where also analytical studies of systems with this decomposition have been done, in this work the trace of 
the extrisic curvature and the determinant of the intrinsic metric are considered dynamical variables and are evolved 
jointly with the rest of the system. 



II. SYSTEM II 

The conformally-decomposed system that we take as starting point has appeared in , and is a variation of the 
system used by Shibata and Nakamura It is a system of 15 equations for 15 variables (</), iiT, 7ij , , f *), and is 
referred to as System II in ||ll|, to distinguish it from the standard 3-1-1 Einstein equations ||T^. These variables are 
related to the intrinsic metric 7^- and extrinsic curvature Kij as follows: 

(la) 
(lb) 
(Ic) 

\l^,K^ (Id) 

(le) 

where 7*-' is the inverse of 7^ . The Einstein equations in terms of these variables are equivalent to the following: 

j^<^ = -\aK (2a) 
^7y = -2ai„- (2b) 
Ax = -fW^D.a + a (a,A'^ + ^K^^) + ^a{p + S) (2c) 





= det(7y)'/' 






K 
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j^A, = e-'^ {-{D.D.af^ + a{R[/ - Sj,^)) + a{KA, - 2AuA',) (2d) 

-27™^^/3''\m ) . (2e) 

Here a is the lapse function, is the shift vector, and r*j, are the connection coefficients of 7ij. The superscript '^^ 
denotes trace-free part, e.g. Rjj^ = Rij — jijR/3. Indices are raised and lowered with 7*-' and its inverse. We use the 
shorthand notation 

where £(} is the Lie derivative along We have, as well, 

i?r = % - p^'i'^Rki (2g) 

i?,, = % + Rf^ (2h) 
Rf^ = -2AA-0 - '2l^Jl'''DkDl<P + iDr^Dj<j) - 4%7'''^fc'/'A</' (2i) 
R^3 = -^7'''7y,fci +7fc(^f^J) 

(2j) 

f = ^7*^' (7.ij + 7jM - 7y,0 • (2k) 

This system is first-order in time and second-order in space, thus it is of second order overall. We show how System 
II can be handled in order to be turned into a well-posed form. Firstly we reduce the system to a straightforward 
first order form, and subsequently we densitize the lapse and combine the constraints into the evolution equations. 



A. System II reduced to first-order form 

We define a set of 12 additional variables 

3 1 

Vijk = lij,k - -7fe(i7j)«,s7"'* + g7^i7fcr^,s7"^ (3a) 

which is trace free in all its indices, namely: Vijkj^-' = and Vijk'y-''^ = Vijk'y^'' — 0, and another set of 3 additional 
variables 

Qi = </>,i • (3b) 

Evolution equations for these new variables are obtained by taking a time derivative of (H) and commuting time 
and spatial derivatives in the resulting right-hand sides. The complete system of equations is now 

j^^=-laK (4a) 

^7y = -2aiy (4b) 

j^K = -^fW,D,a + a (^A,A'^ + ^K'^ + ^a{p + S) (4c) 

= {-{D.D.af' + a(i?,7 - Sf/)) + a{KA,, ~ 2AaA',) (4d) 
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~ 6 ~ 2 - 

,fc "^^»j7 + P' n Vijk + „ l^^ifc - 2a,fe i„- (4f ) 

+ ^7fe(i^")a,n +^Q!^fc(if j) - ^aiijTfe (4g) 
3 1 

+7;(*/3'u)fc -T^P^Ai %k + -flioP\ik (4h) 

-P\ns f^lk(^lj)l ~ ^7*i7fci^ (4i) 

+ (2ai"^ + 7""r ) {^lk{^l,)n,s - ^7y7fc.,n) (4j) 

Q^ - P'Q^,l = -laK,, Qi - ia„ K (4k) 

D D 

where, as before, Rfj^ = Rij — ■|7*-'7'^'-Rfci, with Rij = _Rij + i?^-, and 

i?f^. = -2b,Qj - 2%^'''DkQi + AQ.Qj - A^^'^'QkQi, (5a) 



1 7 ~ 1 ~ 

-Ry" = -T^Vijk, ' +T7T7fe(jr ,j) +T7T7ur ,fe 



f = ^^'(.,) - - ^'jf.f,) + ^7.,f ^ (5c) 

and indices are raised and lowered with 7*-' and 7^^ respectively. The derivatives of the form 7y,fc that appear in the 
right-hand sides of (^) must be interpreted simply as shorthands for combinations of the fields f* and Vijk, via 

3 - 1 ~ 

lzj,k = Vijk + g7fe(»rj) - -lijTk (6) 

For this first-order system to be equivalent to the Einstein equations (in the sense that its set of solutions is the 
same as that of the Einstein equations) , the following set of constraints must be imposed on the initial data (and are 
subsequently preserved by the evolution, as will be shown in the next section): 

n = Y'R^j - A,A'^ + _ 2p (7a) 

P, = ^''DiAj ~ ^D,K + AQiA\ + ^KQ, - S, (7b) 
g- = r^f^„ (7c) 

Q^ = Q^ - (7d) 

3 1 

V.y7c = V^jk - lio.k + g7fc(j7j)«,s7"^ - -7^i7fcn,s7"^ (7e) 



where 



2 ' 10' — 5 

-^Pf + f '^f^^fc + 2f '^""f + f ™^'r„fc, ) . (7g) 
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Constraints (7a) and (7b) are the hamiltonian and momentum constraints of the 3+1 decomposition of the Einstein 
equations, written in our choice of variables. Constraints (|7c|) , (7d) and ( [t^ arise in turning the original second-order 
system into first order. 

In (|)and (|), the derivative Di is the covariant derivative with respect to 7^j, and is related to Di by undifferentiated 
terms: 



rk 
i 



(8) 



B. Taking advantage of the available freedom 

In this section we take advantage of two facts that have been used successfully in similar problems [|l7|,^| j7| . Firstly, 
we densitize the lapse a (and in doing so we introduce a free function, referred as "slicing density" in 

a = e^'"*'G (9) 

Like the shift vector the lapse density u will be considered arbitrary but fixed, a source function independent of 
the dynamical fields. 

Secondly, the evolution equations can be combined with the constraints without altering the set of solutions. We 
add the scalar constraint with a factor h to the evolution equation for K and we add the vector constraint to the 
evolution equations for F* and Q.;, with factors of c and d respectively. In this manner we obtain a system of the form 

u = A\u)\I^u^ B{u). (10) 

A system of this form is known to be well posed if the matrix- valued vector admits a symmetrizer, namely, 

a positive definite, symmetric, bi- linear form ff, in the space of the fields m, whose product with A* (it) yields a 
symmetric-bilinear-form- valued vector. Thus, in order to determine well-posedness, it suffices to consider the principal 
part of the system. 

In this case, the principal terms are 

= (11a) 

= (lib) 

k = P^K.i ~a{Aa + Sh)e"^'l'^''^Qk,i + afoe-^-^f (11c) 



(lid) 



V = l3't\i+acA'\i~'^{c + 2)af'K,i (lie) 
~ 6 - 2 - 

V^Jk = I3%jk,l - 2aA,^^k + ^a7fc(j^j)m,n7"" - ^a7«i^fcm,n7""' (111) 

= /3'Qm - ^(1 + ^c^)-^"- +da^''A,j,i- (llg) 

Our aim is to show that there exist choices of the numerical factors a, b, c, d such that the system ( pl| ) is symmetrizable, 
therefore, well posed. We establish this result by defining a candidate symmetrizer H given as 

p-40 

uHu = (p^+ 5'^5^^%%i + nle~^'*'K^ + A'^ A,^ + n^e^-^-^FT, + -^V^jkV'^'' 

+nle-^'f'Q'Q, (12) 

where ni, n2, are any fixed real numbers different from zero and bounded, so that C^^I < H < CI where C is a 
positive constant and / is the identity operator on the space of u . 

We can easily arrange the values of a,b,c,d so that symmetry of HA^{u) is attained. To this effect, they must 
satisfy 
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rijC (13a) 



-2{2a+l)^nld (13b) 
-(4a + 86)^2 = -1(1 + 4d)nl (13c) 

nlb^~^nl{c + 2) (13d) 

There is clearly plenty of freedom in the choice of a, b, c, d, since any choice that results in non- vanishing ni, 712, 
is allowed. The freedom is thus parametrized by the values of 121,122, n^, since these can take independent positive 
values. Thus our four parameters a,b,c,d are not all independent, but there is a relationship between them that 
reduces the freedom to 3 independent parameters. We can solve (O) for a, b, c, d in terms of 711,122,123, which yields: 



9/5 + 8n^ + n§/8 
(3n2 + 2) 
2/7 



(14a) 



7 



10n2 



(14c) 



It is clear from ( p^ that a and c will take only strictly positive values, and b and d will take only strictly negative 
values, for all real values of 111,122,123 different from zero. 



C. Structure of characteristics 



The system IhV) is of the form 



(15) 



The characteristic covectors are covectors (,a = (Cii ^'^) such that Ci^jl^'' — 1 a-nd such that 

det(A-C,) - (16) 

The values of v that satisfy (|l^) for every direction S,i are the characteristic speeds in that direction. In order to 
find these values we set up an eigenvalue problem for the principal symbol A"^a and find the null eigenvectors. The 
eigenvalue problem is 



n''^a<l> = 
n^^ajij = 

n'^^aK - -(4a + 8b)e-^^^''Qk + foe-^-^^zf 



-2(2a + l)e-4^ (^CoQ,) " ^7y^'^Qfe) 



n°Car'=cCiA*'--(c + 2)fi^ 



^^{i + id)i,K + deAj 



(17a) 
(17b) 
(17c) 



(17d) 

(17e) 
(17f) 
(17g) 



Clearly, n°'£,a = allows for 18 eigenvectors. This is because (17c), (17i) and (17g) in this case constitute an 
overdetermined system of 18 homogeneous equations for 6 unknowns (Aij, K), with zero as the only solution, whereas 
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(17c) and (17d) constit ute a system of 6 equations for 18 variables, which leaves out 12 of the 18 fields (Ky7j,r*,Q) 
free. Lastly, (17a) and (17b) leave the 6 variables (7^ , (f) free. If we represent the eigenvectors in the form 



(18) 



where f(^) := f'^., ff ^ f, - e'^^^.f^^k. V^^^ := V,,ke,vg^ - e-^^^kV^.i^, ^^^^^ A""^^ ■= 

A^j£_3 - e-^'t'^^A''^^k^l, A^J'^^ := A^ - 2e-'*'^C(iij)/C' + e'^'^iiiiA^^ikiu then we have 5 eigenvectors corresponding 
to the five components of the conformal metric 

(%-, 0,0, 0,0, 0,0, 0,0,0); 

we have the determinant as an eigenfield 

(0,0,0,0,0,0,0,0,0,0); 
we have 7 eigenvectors corresponding to the seven transverse components of Vijk 

(0,0,0,0,0,0,y,g\ 0,0,0); 
we have 3 eigenvectors corresponding essentially to the three components of Qi 



(0, 0, g„ ^^^rQ„ 0, ~2(2a + l)e(.Q,) - i% (^^^^^^ - 2(2a + 1)^ Q^^', 0, 0, 0, 0); 

and we have 2 eigenvectors corresponding essentially to the two components of the transverse part of Vi 

(0,0,0,0,ff\^e(.f5p,0,0,0,0); 



(19a) 
(19b) 

(19c) 
(19d) 

(19e) 



If n°'^a 7^ 0, then 7^ =0 = 0, and we can solve (17e), (171) and (17g) for (V^y/c, F', Q) in terms of and (A^ , K). We 
can substitute (V^ij^, f Q) into (17c) and ( |l7d| ), obtaining thus a system of 6 equations for the 6 variables {Aij,K) 
as follows 



= Ke^'^ (^{n^Uf - i(4a + 86)(1 + 4d) + h{c + 2) j + ^-A-^ ((4a + 8b)d - be) 
= A,e^^ (1 - (n^Caf) ~j(^l{c + 2)- (2a + 1)(1 + 4d)^ (^^.0 - ^e^^^ 



,10 5 

where we have used the notation 



^ - ^ - 2d(2a + 1) ) ( r - -7yf 



(20a) 



(20b) 



(21) 



If 1 — (ri^^a)^ — 0, then K = (,-A^ = by (20a), which implies ^'^h = by (20b). However, two of the five components 
of Aij are thus free, which means that there are 4 eigenvectors, essentially labeled by the transverse components of 
Aij. We have 2 eigenvectors for n^^a = 1 



and 2 eigenvectors for n°'^a = ~1 



(o,o,o,o,o,-2l^^\ 0,0, 0,0,1^^0; 



(0,0,0,0,0,2i|J-\0,0,0,0,A^J-) 



(22a) 



(22b) 



If 1 - (n^^a) ^ 0, then contracting ( |20b| ) with C yields 
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^ ^!£-^-2rf(2a+l)) -KCaf 



-le^'^K^j Q(c + 2)-(2a + l)(l + 4rf)^ 



1 



6 



3 



+ ^&A-a, -c---2di2a+l) 
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(23) 



Thus, if H-i(7c/10-3/5-2d(2a+l))-(n"^a)^ = 0, then K = ^-A^ = by ( pOa| ), which imphes that Q is identically 
satisfied, thus three out of the five equations (20b) are identities, the remaining two determining two components of 
Aij. Thus two of the five components of Aij are free, which means that there are 4 eigenvectors, essentially labeled by 
the two longitudinal-transverse components of Aij . We have 2 eigenvectors for n°'^a = \/ (3/5 — 7c/10 — 2d(2a + l))/2, 
namely 



d ~ 



(0, 0, 0, TT^r^ ^ ^ ( %(.^ 



C 7(LT) 



(LT) 



6 



Ci ' ' ' Ci 
and 2 eigenvectors for n°'^a 



5Ci 



5Ci 



J(LT) 



(LT) 



0,0,if^\0), (24a) 



V(3/5 - 7c/10 - 2d(2a + l))/2, namely 



(0,0, ^^A, ,0, ^^A, ,^^^^i^A^) 



(LT) 



6 
5CT 



-.(LT) 



(25) 



where we have used the shorthand notation 

Ci V(3/5 - 7c/lG - 2d{2a + l))/2 

But if 1 + i(7c/lG - 3/5 - 2d(2a + 1)) - {n''£_a)^ ^ 0, then ^^Aij is determined by the values of K and ^-M by (H), 
and if plugged back into (20b) it follows that all the components of Aij are determined by K and ■C'^"?- Therefore it 
is necessary that K and ^-A-^ be nonvanishing. Contracting ( |23| ) with we obtain 

2 







^e^^i^ f^(c + 2)-(2a+l)(l + 4d) 
) \5 



H-A-ai 



7 

3 Vio' 



2(2a+l)d - (n%) 



(26) 



Equations (P0a| ) and (26) form a system of two homogeneous equations for K and ^-A-^. Thus, for K and ^-A-^ to 
be nonvanishing, it is necessary that the determinant of the system be zero. The determinant is 

— (-3(n°Ca)^ + 2a)(15(n"^a)^ - 9 + 105c - 80bd + 2Qd - 7c) (27) 
45 

It can be seen that, because a and c are strictly positive and b and d are strictly negative, the four roots of the 
determinant are real. For the roots n°'^a of the determinant, we have 



K = - 



91_ 
2 



2 (7c 

3 V 10 



2d(2a+ 1)) 



|(c + 2) - (2a + l)(l + 4c;)) 



-40 



Q^ = 



_ 3(1 + M) 1 - {n-^a f + f - f - 2d(2a + 1)) g-^^ 



|(c + 2)-(2a+l)(l-F4d)) 



U-A-^ 



c-3{c+2 



1 - {n-^a? + 1(^-1- 2d(2a + 1)) \ ^-A-^ 



\{c + 2)~{2a+l){l + 4d)) 



r 



(T) 

i 

(L) 
(T) 



5n"^a 



ijk 



= 



-40 



(28a) 
(28b) 

(28c) 

(28d) 
(28e) 

(28f) 
(28g) 
(28h) 
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These are clearly four distinct eigenvectors, since there are four distinct values of n"^a given by the roots of the 
determinant (^). These can be thought as being labeled, essentially, by K or £,-A-^ indistinctly, or we can associate 
one characteristic speed to K and the other one to ^-A-^, as we prefer to do below. 

Summarizing, we have characteristic speeds obtained from the following distinct values of n'^S,a' 

^) n°-^a — 0, timchkc, with eigenfields (essentially) jij , (j), Qi,f^'^\V^T' . 



t>) 'n°'£,a = 1, null, with eigenfields (essentially) Jl^"^"^-* 



(TL) 



c) n^^a = (1 + 5(7c/10 - 3/5 - 2d(2a + 1)))^/^ = Ci, with eigenfields (essentially) A 

d) n"^a = (2a/3)i/2 ^ (j^^ ^^^-^ eigenfield (essentially) K 

e) n"Ca = (3/5 - 2foc/3 + 165(i/3 - 4d/3 + Tc/lS)^/^ = C3, with eigenfield (essentially) i^^^) 

In the expressions for the characteristic speeds c),d) and e), the parameters a, 6, c, d are given in terms of ni, 712, na 
via (|l4[). These speeds may be superluminal or causal, depending on the values of ni, rt2, ^3- We can choose ni, 712, na 
so that Ci , C2 and C3 are all equal to 1 . This is achieved by setting 

^2 _ 4 280n2 + 49 + 400n| ^^9^) 



^ 15 mnl - 49 
2 6400ni 

for any value of n| greater than 49/60. This means that there is a one-parameter family of well-posed conformally- 
decomposed systems with "physical" characteristics. From the analytical point of view, there does not appear to 
exist an argument for singling out a preferred value of n2- It is likely that a preferred value of rt2 will be dictated by 
optimal numerical behavior. The expressions for a, 5, c, d in terms of n2, with ni and na as above (|2^), are as follows: 

3 

a = -, (30a) 

h^-''f^-% (30b) 
4(7 + 20^2)' ^ ' 



60n2 - 49 

d^ (30d) 

800nl ^ ' 



D. Propagation of the constraints 

The propagation of the constraints can be calculated by taking a time derivative of each one of the constraint 
expressions, and subsequently using the evolution equations ( [ll| ) to eliminate the time derivative of the fields in the 
right-hand side in favor of spatial derivatives, which recombine to yield back the constraints. We obtain 



n = (3^H,i +(c - M)ae~^'*'')^^Vkj + ■■■ (31a) 

-2a(2a+l)e-^V'Q[Mb + --- (31b) 

g'^p'g\i + --- (31c) 

Q^=^3'Q^,l + ■■■ (31d) 

V,,k = P%k,i + ■■■ (31e) 



where • ■ • denotes undifferentiated terms proportional to the constraints themselves. To analyze the constraint prop- 
agation we proceed to turn ( |3l| ) into first order by defining several sets of variables which represent all the spatial 
derivatives of Vijk,G^ and Q^: 
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— V -, ^ 




^ijkl 




1, 






















= ^{1,3) + 



(32a) 

(32b) 

(32c) 
(32d) 
(32e) 

(32f) 

where A is a constant which will be fixed shortly. Calculating the time derivative of these we obtain the resulting 
first-order system of evolution of the constraints: 

ii = (3^n,i +aic - M)e-'^'l'Vi! + ■ • • (33a) 

A = |3'V^.l + ^(1 - 46)H„ -^e-^^Wa! -2a(2a + 1)6-4"^^,;/ 
2 

~Aae-^^Zu! +^^"4^ J,/ + ^-^e~'^Tu! + • • • (33b) 
mj = - + + • • • (33c) 

■^yfci = P™XijU,m + • ■ • (33d) 

Z^y = + j] + • • • (33e) 
Z,j = + + • • ■ (33f) 

T^j = 13% + acP[,,,] + • ■ • (33g) 

j^j = /3' + a (^c + yAd^ + • • • (33h) 

5' = P'Q\i + • • • (33i) 
= /5'Szj + --- (33j) 
%k = l3'V^jk,i + ■■■ (33k) 

For this system there is a symmetrizer given by 

uHcu = e^^ ^j-^^l n' + V.V^ + e-^^lw., W^^" + X^jkiX^^'^' - 2(2« + 1) 
6(c — 8a) 2c d 

rf^^^^ +^ 30c ^^^^ +^ c + Ad30/7^'^^ 
+Q^Q' + Q^Q' + V^,kV''''. (34) 

Taking A = ^(— c/d), iJc is positive definite because, under the conditions ([ij), all the factors accompanying the 
squares of the fields are strictly positive. This shows that no additional restrictions on the ranges of the parameters 
a, 6, c, d are necessary in order to have well posed constraint evolution. 



III. CONCLUSION 



We have derived a 3-parameter family of well posed versions of the conformally-decomposed 3+1 equations, perhaps 
amenable to successful numerical integration. One might object that there is no need for it in view of the results in pl| , 
but we can argue rather strongly that these results may prove helpful in pinning-down the main cause of numerical 
instabilities. This well posed version requires the lapse to be proportional to the determinant of the intrinsic geometry 
of the surfaces, and requires combinations of the contraints with the evolution equations. The lapse density a and 
the shift vector are arbitrary non-dynamical variables, which means that they must be specified as free source 
functions. This well posed version uses the same variables as the original system (except for the addition of the first 
spatial derivatives of the densitized 3-metric, referred to as "conformal metric" by the authors in ||l^). In addition, 
this well-posed version of the original equations propagates the constraints in a stable manner, which is relevant to 
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unconstrained evolution. We think that this is the least invasive way to turn the original conformally-decomposed 
system into a well posed one. In practice, a choice of the numerical parameters ni,n2,n3 must be made. The 
characteristic speeds depend on this choice. 

Optimal choices of the parameters 711,712,^3 for numerical evolution arc those that ensure that the characteristics 
are all either null or timclike. With such a choice, the formulation would be suited to evolve blackhole spacetimes 
outside the event horizon. Among these choices, it has been suggested that the preferred one would be the one 
for which the characteristics are all "physic al", na mely, either null or normal to the slices. We have shown that such 
a choice is possible for an arbitrary 712 > 1/49/60. 

The systems obtained in this work are not contained in our previous work The choice of variables in 

is inadequate for decomposing the trace and trace free part of the extrinsic curvature, as well as for extracting the 
determinant of the 3- metric. This is clear from the fact that, in that work, the available parameters a and /? are not 
allowed to take the value —1/3 without the argument breaking down. 

The systems obtained here differ significantly from the system obtained in |l^ by considering the trace of the 
extrisic curvature K and the determinant of the instrinsic metric (and its derivatives) as dynamical variables on equal 
footing with the rest, rather than as free source functions. Furthermore, we have obtained a 3-parameter family of 
systems, one system for each appropriate choice of ni, n2, 713, whereas in p5| there is only one system which preserves 
the trace conditions. Additionally, we have separated the divergence of the intrinsic metric F* from the divergence- free 



part of the metric. This decomposition keeps up with the spirit of |11|. 

We have found that in obtaining these well-posed formulations the lapse must be proportional to some power of 
the determinant of the intrinsic metric, since the parameter a cannot take the value 0. This is similar to our findings 
in as well as other notable cases |^,^,|l],0. In our present case this is remarkable, since we have used quite 

general energy norms. 
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